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CHARACTERISTICS OF AEROELASTIC INSTABILITIES IN TURBOMACHINERY - 
NASA FULL SCALE ENGINE TEST RESULTS 


Joseph F. Lubomski 
NASA-Lewis Research Center 
Cleveland, Ohio 


Abstract 

Several aeromechanical programs have been con- 
ducted in the NASA./USAF Joint Engine System Re- 
search Programs. The scope of these programs, the 
instrumentation, data acquisition and reduction, 
and the test results are discussed. Data pertinent 
to four different instabilities were acquired; two 
types of stall flutter, choke flutter and a system 
mode instability. The data indicates that each in- 
stability has its own unique characteristics. 

These characteristics are described. 


Nomenclature 

i 

H 

C-| negative damping coefficient 

C 2 positive damping coefficient 

G shear modulus of elasticity 

I 0 mass 'moment of inertia 
J polar area moment of inertia 
e constant in Van der Pol's equation 

Z damping coefficient 

6 dimensionless angular displacement 

6 angular displacement 

oif frequency of unstable motion 

natural frequency of blade or vane 

Introduction 

Problems due to aeroelastic instabilities 
have occurred in recent development engines that 
have resulted in costly and time-consuming program 
delays. These have provided an impetus throughout 
industry and government to develop a better under- 
standing of the nature of aeroelastic instabilities 
and to evolve design tools to avoid such problems. 
The NA.SA. and the USAF have several joint programs 
in this effort. One of these is the Joint Engine 
System Research Programs, the objective of which 
is to investigate technology areas of mutual inter- 
est, one of which is aeroelastic instabilities. 

In the Joint Engine System Research Programs 
a thoroughly instrumented engine is installed in 
an altitude facility and tested at simulated 
flight conditions. For the aeromechanical programs 
the flight conditions selected were those at which 
aeroelastic instabilities were known to exist 
based on engine development histories. Instabil- 
ities were induced at these test conditions and a 
comprehensive set of aerodynamic and aeromechanical 
data were obtained prior to entering and during 
the instability regions . This data is now a part 
of a data bank available to help understand aero- 
elastic phenomena and for use in developing and 
verifying design, tools. This data has been vali- 
dated and does define the behavior of aeroelastic 
instabilities in an engine at flight conditions. 

To date two sets of aeromechanical tests have 
been completed. Data relative to stall flutter, 
choke flutter and a system mode instability has 
been collected. This was accomplished using two 


different engines, an afterburning turbofan engine 
and an afterburning turbo-jet engine. Both en- 
gines are high thrust to weight ratio, high per- 
formance engines representative of current tech- 
nology. 

This paper will outline the scope of these 
two test programs, the test procedures and in- 
strumentation used, the data acquisition and re- 
duction methods used, and the test results. 

Types of Aeroelastic Instabilities 

An aeroelastic instability is a self-excited 
vibration of an airfoil, single or cascaded, 
whereby the unstable motion is generated and sus- 
tained due to the unique motion of the airfoil in 
the airstream. Thus the oscillating airfoil is 
extracting energy from the airstream and is aero- 
dynamically negatively damped. 

Aeroelastic instabilities in turbomachinery 
are generally now classified as being one of five 
distinct types-*- relative to the location of the 
instability on a compressor map such as shown in 
Fig. 1. These are as follow: 

I. Subsonic/ transonic stall flutter. 

II. Choke flutter. 

III. Low back pressure supersonic flutter. 

IV. High back pressure supersonic flutter. 

V. Supersonic stall flutter. 

All of these instabilities are being investi- 
gated in current NASA/USAF programs. However, in 
the Joint Engine System Research Programs to date, 
no supersonic instabilities have been encountered, 
thus the discussions to follow will be limited to 
subsonic/transonic stall flutter and choke flutter 
with one exception. This exception is a system 
mode instability relative to which only a limited 
amount of data has been acquired thus making it 
inappropriate at this time to attempt to classify 
it unequivocally into one of these types. 

Subsonic/transonic stall flutter was histor- 
ically the most commonly found aeroelastic in- 
stability in turbomachinery prior to the develop- 
ment of high tip speed fans operating supersoni- 
cally at the tip. It occurs near the stall line 
at large angles of incidence and is accompanied 
by flow separation for part of or all of the os- 
cillatory cycle. The flutter frequency is near 
but not necessarily equal to one of the airfoil 
natural frequencies in one of the first three 
fundamental modes, viz; first bending, first tor- 
sion or second bending. 

Choke flutter occurs near the choke line at 
small positive or negative incidence angles and is 
accompanied by choking in the cascaded blade pas- 
sage. The flutter frequency is again near but not 
necessarily equal to a blade natural frequency, 
usually first bending. In earlier engines choke 
flutter usually affected turbine stages rather 
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than compressor stages and is still a problem on 
high aspect, low pressure steam turbines. 

A. system mode instability^ will normally 
afflict a shrouded blade rotor assembly. Insta- 
bilities of this type are dependent on shroud 
coupling and occasionally disk coupling. They are 
characterized by flutter frequencies not necessar- 
ily associated with a fundamental blade frequency 
and the blade vibratory made shape is more often 
than not a coupled mode, such as a coupling of 
first bending and first torsion. This type of in- 
stability has not been thoroughly documented. 

Apparatus 

The two series of tests completed to date 
were conducted using two different types of en- 
gines. The first series of tests were done using 
a twin spool, high performance, afterburning turbo- 
fan engine which had a history of an aeroelastic 
instability with a preproduction first stage fan. 
The installation of this engine in the test facil- 
ity is shown in Fig. 2. This engine has a three 
stage fan as shown in Fig. 3. The fan blades have 
part span shrouds and are made of titanium. The 
bypass airflow ratio is nominally 0.7 to 1 at 
sea level static standard day conditions. The 
diameter at the fan inlet is nominally 88.34 cm 
(34.78 in.). This engine was susceptible to 
subsonic/transonic stall flutter in above shroud 
torsion at elevated inlet temperatures and pres- 
sures i.e. at high Mach number flight conditions. 
This particular engine was tested in the PSL-1 al- 
titude chamber of the Lewis Propulsion System Lab- 
oratory. 

The second engine tested was a single spool 
afterburning turbojet. The installation of this 
engine in the test facility is shown in Fig. 4. 

This engine has a 9 stage titanium bladed com- 
pressor as shown in Fig. 5. The first two stages 
have part span shrouds. The inlet guide vanes and 
the first three stators are variable. The inlet 
diameter of this engine is 43.61 cm (17.17 in.). 
This engine was tested in the PSL-2 altitude cham- 
ber of the Lewis Propulsion Systems Laboratory. 
This engine had development problems with several 
types of instabilities. The third stage and the 
fifth stage rotors were susceptible to choke flut- 
ter. The fifth stage rotor instabilities occurred 
at cold inlet temperatures, less than -20° C 
(-4° F) . The third rotor instabilities occurred 
over a range of inlet temperatures bounded by 
15.5° C (60° F) and 49° C (120° F) at elevated in- 
let pressures. The second rotor had exhibited a 
system mode instability at elevated inlet pres- 
sures and temperatures. The fourth stage stator 
had a stall flutter instability at high inlet 
pressures and cold inlet temperatures. 

Instrumentation 

Data in the aeromechanical programs were ac- 
quired using a wide variety of instrumentation and 
data reduction techniques. Jones, et al.,^ com- 
prehensively discusses the type of instrumentation 
and the location of this instrumentation in the 
engine as used in the turbofan program. The 
turbojet program had essentially the same types of 
instrumentation and used the same data acquisition 


and reduction techniques. As background informa- 
tion the aerodynamic and aeromechanical instrumen- 
tation will be summarized herein. 


Aerodynamic Instrumentation 

The principal measurements taken are total 
temperatures, total pressures and static pressures. 
These are obtained at all principal engine sta- 
tions using a number of radial rakes. Wall static 
pressures are also obtained at these stations. 
Interstage data was also obtained using 5 element 
radial temperature and total pressure rakes along 
with two wall static pressure measurements at each 
location. In the turbofan program this instru- 
mentation was mounted into the leading edge of a 
stator vane in each of the three fan stator stages 
to avoid blockage. In the turbojet program probes 
were designed and built that fitted in the stator 
rows for each of the nine stages. A number of 
these were strain gaged and monitored during the 
test program. The arrangement of the aerodynamic 
instrumentation by stations for the turbofan pro- 
gram is shown in Fig. 6. The arrangement in the 
turbojet program was similar. 

The data obtained from this instrumentation 
was recorded on two systems. One system is used 
to record the steady-state data points and the 
other system is used to record transient data. In 
the transient system data is obtained during a 
transient with a sample rate of up to 150 a second. 
This recorded data is then processed and reduced 
off-line using a large scientific computer. A 
significant accomplishment of the aeromechanical 
programs was the acquisition of comprehensive 
aerodynamic data while an instability was occurr- 
ing in the engine. Normally the instability was 
sustained at a predetermined safe stress level for 
from 1 to 2 minutes while data was being acquired 
on all data systems . 

Aeromechanical Instrumentation 

In the aeromechanical programs data was ac- 
quired using three types of instrumentation. The 
primary source of data was strain gages mounted 
on the fan and compressor blades and vanes. The 
data from these were used for on-line monitoring 
as well as for post-run analysis. In addition a 
photoelectric scanning (PES) system was de- 
veloped^ which used light probes. This system 
made it possible to observe online the nonintegral 
response of all the blades on a given rotor. The 
-data from this system was also recorded as -either 
video tape or digitized for post run analysis. A 
third method, developed by Kurkov^ utilized data 
recorded from high response pressure transducers 
which when analyzed in the frequency domain proved 
to be useful in defining the aeromechanical modes 
of a rotor assembly in flutter. 

Strain gages were mounted on selected blades 
and vanes by the engine manufacturers after bench 
testing the blades and vanes to determine the 
exact natural frequencies and mode shapes . In the 
turbofan program all of the fan blades were bench 
tested individually and as an assembly to deter- 
mine the system modes . In the turboj et program 
only a select number of blades were bench tested 
due to the large number of blades and vanes in- 
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volved in the program. In the turbofan program the 
three fan rotor and stator stages were strain gaged 
with the bulk of the strain gages being on the 
first stage rotor since only it was susceptible to 
an aeroelastic instability. Three fan builds were 
made and tested and in each build there was a 
slight variation in the strain gage distribution. 
Only one build was made and tested in the turbojet 
program. In this build all nine stages had strain 
gaged rotor blades as did the last six stator 
stages . A 100 channel slipring was used in both 
programs to transmit the rotor strain gage signals 
thus 50 rotor strain gage signals were monitored 
at a time. In addition to these between 30 and 40 
stator gages were monitored and recorded. 

Strain gage signals are monitored in the time 
domain and in the frequency domain. Each strain 
gage signal is displayed on a small oscilloscope 
for continuous monitoring in the time domain. 
Currently up to 112 signals can be handled simul- 
taneously. These signals can also be quickly 
patched into large oscilloscopes for expanded 
analysis in the time domain or they can be se- 
lectively patched into a number of systems for 
analysis in the frequency domain. In the aero- 
mechanical programs the following systems were 
available and were used: 

1. A software based Fast Fourier Transform 
(FFT) Digital Signal Analyzer which can handle two 
channels simultaneously. 

2. Several analog spectral analyzers. 

3. A microprocessor based multichannel vibra- 
tion and flutter monitor® which can handle up to 
16 channels simultaneously. 

All strain gage data is recorded for post run 
analysis. Forty FM channels are available within 
the test facility for the more pertinent strain 
gages and 180 FM multiplexed channels are available 
in a central analog facility for all the strain 
gages. Reduction of strain gage data is done in 
the time domain and in the frequency domain. Time 
domain data reduction is done to obtain overall 
stress levels, waveform and event history. Fre- 
quency domain data reduction is used extensively 
due to the availability of several FFT systems. 

In the frequency domain the data is processed to 
obtain power spectral density, linear spectra, 
transfer functions and cross correlation function. 
The coherence function is always computed whenever 
time averaging is used. 

Results and Discussion 
Turbofan Program-Subsonic Stall Flutter 

The program using the turbofan engine was con- 
ducted in three phases during which 33 data points 
in flutter were obtained at a variety of inlet 
conditions and speeds. Aerodynamic and aeromech— 
anical flutter data were obtained at inlet pres- 
sures between 103 and 207 kfa (15 and 30 psia) and 
at inlet temperatures from 135° to 171° C (275° to 
340° F) . The flutter data was obtained at cor- 
rected fan speeds between 64 and 73 percent of 
design speed. The first phase of the program was 
conducted to develop monitoring techniques, in- 
strumentation and operational procedures and to 
locate the flutter regions within the engine 
flight envelope. The instability region was 
approached by setting inlet conditions and speed. 


Then the nozzle was throttled while holding con- 
stant corrected speed to increase fan back pres- 
sure. Only stall flutter in above shroud torsion 
of the first stage fan rotor was found and docu- 
mented at several conditions. The fan assembly 
was then removed, disassembled, instrumentation 
refurbished and the fan assembly rebuilt for 
Phase II testing. It was during the Phase II test 
program that the majority of the data was obtained. 
For the Phase III build the instrumented fan blades 
were removed and replaced with different instru- 
mented blades. Using this build several of the 
previous test points were repeated and additional 
flutter data was obtained which proved to be of 
significance. 

The response of the rotor blades in flutter 
was varied. At incipient flutter a low level ir- 
regular beating was observed on the strain gage 
signals probably due to flow separation. Upon 
developing full flutter the strain gage signals 
usually became periodic. Several flutter fre- 
quencies were observed of which one was predominant 
on most of the strain gages. A blade to blade 
variation in maximum amplitude was observed on the 
strain gages and on the light probes. Usually a 
small group of blades would respond initially upon 
penetrating the flutter boundary and then as the 
excursion into the flutter region progressed more 
of the blades would become active. 

No unique and common interblade phase angle 
was observed from the strain gage data. Light 
probe observations indicated that a few small 
groups of blades may. have been coupling with a 
common interblade phase angle but complete coupling 
of all the blades on the rotor assembly did not 
occur. The excursion into the flutter region was 
normally stopped when any single strain gage in- 
dicated an overall stress level of 75 percent of 
allowable, as specified by the engine manufacturer. 
A few excursions were made beyond this to see if 
complete coupling could be induced but were un- 
successful. Thus the flutter or aeromechanical 
mode appears to be mistumed as defined by White- 
head . 7 

The flutter frequency observed usually did not 
vary while the rotor blades were in flutter as 
shown in the expanded spectrum plot of Fig. 7. 

This data is from a strain gage located slightly 
above the shroud at the trailing edge. This blade 
is responding at 1072 Hz, which was the primary 
flutter frequency, with a maximum amplitude of 
±50.3 MPa (±7.3 ksi) at that frequency. The re- 
sponse at 1040 Hz is a secondary flutter fre- 
quency from a blade fluttering at that frequency 
located three blades away on the rotor. This will 
be discussed in detail later. The response at 
1120 Hz is the eighth engine order response. Dur- 
ing the time that the rotor was in flutter these 
frequencies did not change and the interblade 
phase angle between any pair of adjacent blades 
remained highly coherent. 

An exception to this behavior was observed at 
either the lower inlet pressures or inlet temper- 
atures whenever flutter occurred. At these inlet 
conditions the flutter frequencies varied irreg- 
ularly as did the amplitudes. This is shown in 
the expanded spectrum plots of Figs. 8 and 9. 

This data is from the same strain gage as dis- 
cussed previously but for a lower inlet temper- 
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ature, 135° C (275° F) . These data are single 
average FFT linear spectra taken 1 minute and 20 
seconds apart during the same steady-state data 
point. In between this time there were no signifi- 
cant changes in engine parameters. These data 
show a number of flutter frequencies occurring one 
of which is common to both time frames at 1064 Hz. 

A frequency content such as those shown will yield 
an erratic beat response in the time domain com- 
monly associated with flow separation vibrations. 

A, deeper excursion into the flutter zone was made 
once to attempt to develop a unique flutter fre- 
quency but this only resulted in rotating stall 
which annihilated the aeroelastic instability. 

Time averaging this type of data in the frequency 
domain usually shows between 3 and 4 predominant 
flutter frequencies with good coherence of most of 
the interblade phase angles. A similar response 
was observed in data taken at an inlet pressure of 
103 kPa (15 psia) and an inlet temprature of 
171° C (340° F) . Time averaging of these data 
yielded a poor coherence function relative to the 
interblade phase angles, less than 0.4, and only 
one predominant frequency. Thus it would appear 
that there exists a domain which is a function of 
stage inlet conditions where the flutter response 
is irregular and aperiodic that separates the 
well defined flutter domain and the flutter clear 
domain. 


Inlet pressure and temperature besides affect- 
ing the rotor blade's response also had a signifi- 
cant influence on the flutter boundary location on 
the fan performance map. Over the range of inlet 
pressures and temperatures investigated in the test 
program it was found that increasing inlet temper- 
ature and increasing inlet pressure were both de- 
stabilizing. That is, increasing either of these 
parameters tended to move the instability boundary 
nearer the operating line. This is shown in Fig. 

10 which is a segment of a fan performance map. 

Five data points obtained at approximately the 
same corrected speed from a single strain gage at 
nominally the same stress level are plotted show- 
ing the effect of inlet pressure and temperature 
on the flutter boundary. 

Another aspect of varying the inlet pressure 
and temperature was the effect on the flutter fre- 
quency. Increasing the inlet pressure had the 
most significant effect. This is shown in Fig. 11 
which is a plot of the flutter frequency obtained 
from a single strain gage mounted on the blade in 
position 1 on the fan rotor. Two sets of fre- 
quency data are shown. One set encompasses the 
frequency range between 1060 Hz and 1900 Hz which 
was the primary or most common flutter frequency 
observed on ell the instrumented blades. The sec- 
ond set of data between 1020 Hz and 1050 Hz is a 
secondary flutter frequency which commonly showed 
up only on those blades around position 1 on the 
rotor. It should be noted that the behavior of 
the two sets of data with respect to inlet temper- 
ature and pressure is similar, thus the two sets 
of data are probably a common phenomena. It was 
found that increasing inlet pressure tended to 
clearly decrease the flutter frequency in what 
appears to be a linear function. Also increasing 
the inlet temperature slightly decreased the flut- 
ter frequency. This latter effect Is probably due 
to the temperature effect on the shear modulus. 

Another effect not shown but which affects 


this data was the effect of physical or mechanical 
speed. Increasing the engine speed for a given' 
inlet temperature and pressure tended to slightly 
decrease the flutter frequency. As an example, 
two data points taken at inlet conditions of 138 
kPa (20 psia), 171° C (340° F) and 700 rpm apart 
had a difference in flutter frequency of 4 Hz 
with the higher speed having the lower frequency. 
This was typical of all the data taken. 

The functional behavior of the flutter fre- 
quency with an engine parameter may be highly 
significant. Aeroelastic instabilities can be 
considered to be a limit cycle phenomena. Den 
Hartog® developed a simple single degree of free- 
dom model which can be. analogous to the behavior 
of a blade in flutter. Den Hartog's model formu- 
lated as a single degree of torsion model yields 
the equation: 

I 6 - (C 1 - C 2 6 2 )0 + GJ0 = 0 (1) 


where Cq is the negative part of the damping co- 
efficient and C£ the positive part. As can be 
seen from Eq. (1) small values of 9 will yield 
a negative damping coefficient while large values 
of 9 will yield a positive damping coefficient. 
For periodic motion the sign and value of the 
damping coefficient will vary during a single os- 
cillation. Proceeding with Den Hartog's develop- 
ment and making the variables dimensionless yields 

0 - e(l - 0 2 )0 +0 = 0 (2) 


where 


- w 


(3) 


Eq. (2) is a common form of Van der Pol's equation 
which defines the behavior of a nonlinear, self- 
excited relaxation oscillation. The period for a 
relaxation oscillation is 2e/a n or 



(4) 


Thus as e increases with m n remaining constant 
the frequency of the relaxation oscillation will 
decrease. The Van der Pol coefficient e will 
increase, holding I Q and u n constant only if 
the negative part of the damping coefficient, Cq, 
increases . 


Applying this development to the flutter fre- 
quency data it can be hypothesized that the change 
in flutter frequency with increasing inlet pres- 
sure may be due to the effect of the stage inlet 
pressure on the negative work done during an 
oscillatory cycle. The small change in flutter 
frequency due to the inlet temperature is prob- 
ably due to the change in natural frequency ((%) 
because of change in structural stiffness as a 
result of temperature. In a similar manner in- 
creasing rotor speed tends to decrease the tor- 
sional natural frequency, thus reducing the flut- 
ter frequency, which is consistent with data ob- 
servations. 


The effect of changing the Instrumented 
blades and repeating several data points yielded 
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several interesting results. The secondary flutter 
frequency (1020 to 1050 Hz) did not appear. The 
location of the maximum response on the rotor as 
observed with light probes shifted by about 100°. 
That is the aeromechanical mode of the rotor 
changed, but still remained mistuned. In conjunc- 
tion with this the flutter boundary on the fan per- 
formance map, as established by the strain gage 
on the blade in position 1, moved closer to the 
operating line. That is, in Phase III for the same 
conditions flutter initiated at a lower pressure 
ratio than that observed in Phase II. The pres- 
sure ratio at which flutter was found in Phase III 
was approximately 3 percent lower than in Phase II. 
The flutter frequencies were essentially the same 
varying only by about 3 Hz. 

The disappearance of the secondary flutter 
frequency may be a result of changing an instru- 
mented blade near position 1. It should be noted 
that none of the strain gages on this blade were 
functional making it impossible to verify the 
supposition being made. This blade had a room 
temperature bench test natural frequency of 1045 
Hz, the lowest of all the blades in the rotor 
assembly and thus making it a likely candidate to 
be susceptible to stall flutter. It was replaced 
by a blade with a much higher torsional natural 
frequency. The change in aeromechanical mode and 
the apparent decrease in stall flutter margin can- 
not be readily explained. The most responsive 
blades in Phase II, and the most responsive blade 
in Phase III, were common to both phases. It may 
be conjectured that in making the rotor build for 
Phase III build up tolerances may have changed 
sufficiently to make different blades more sus- 
ceptible to flutter as for example an incidence 
change. 

The concept of establishing a flutter bound- 
ary on a fan performance map using either a single 
or a few strain gages is an arbitrary procedure. 
Mistuned flutter modes are common in an engine and 
one does not know that the instrumented blade be- 
ing used is the most responsive in any one flutter 
mode. Thus the changes in flutter boundary noted 
herein should not be judged quantitatively. 

Turbo j et Program-Choke Flutter 

The majority of the data obtained in the tur- 
bojet aeromechanical program involved the choke 
flutter phenomena. Thirty one data points at a 
variety of inlet conditions and speeds were ob- 
tained documenting this phenomena as it affected 
the third stage rotor. This instability occurred 
in the first bending mode. Comprehensive aero- 
dynamic and aeromechanical flutter data were ob- 
tained at inlet pressures between 103 and 200 kPa 
(15 and 29 psia) with inlet temperatures between 
15.5° C and 49° C (60° and 120° F) . The insta- 
bility occurred at corrected speeds between 95 and 
103 percent of design speed. In order to induce 
the instability it was necessary to off-schedule 
the compressor variable geometry since the oper- 
ational configuration of the engine was free of 
this instability. Flutter data was obtained with 
the first stator opened off-schedule 4.4° and also 
with the second stator closed off-schedule 2° and 
4.3°. Only one stator was off-schedule at any one 
time. Opening stator 1 allows an increase in 
mass flow which tends to slightly decrease the in- 


cidence on rotor 3 while closing stator 2 directly 
decreases the incidence on rotor 3. 

The response of the rotor blades to the in- 
stability was quite different from the stall flut- 
ter response discussed previously. All of the 
strain gaged blades responsed at the same flutter 
frequency but with varying blade to blade maximum 
amplitudes. The instability was approached by 
doing a slow acceleration up along the operating 
line and flutter came in softly. That is, the 
stress levels built up gradually and were con- 
trolled. There was no warning of incipient flutter 
as there is with stall flutter as characterized 
by the erratic beat phenomena. There was a regular 
steady beat more often than not as' the instability 
developed which persisted as the instability grew. 
This was due to the presence -of an engine order 
response in close pro.ximity to the flutter fre- 
quency as seen in the spectrum plot of Fig. 12. 

This spectrum is typical of all the data obtained. 
Common to all the data analyzed was the coexist- 
ence of a blade engine order resonant response 
with the flutter response. The instability fre- 
quency varied with speed, inlet pressure and inlet 
temperature but at no time did it vary more than 
5 percent from an engine order. At times it was 
higher, at times lower, and occasionally was co- 
incident with an engine order. This is quite un- 
like stall flutter. A stall flutter response has 
been observed to occur at an engine order fre- 
quency but more often than not the flutter fre- 
quency will be well removed from an engine order 
response frequency and will be independent of it. 

The effect of increasing either inlet temper- 
ature or inlet pressure was found to be destabi- 
lizing. If either one or both of these were in- 
creased the instability would occur at a lower 
corrected speed. This can be seen in Fig. 13 
which is a plot of flutter data points as a func- 
tion of inlet temperature versus mechanical speed 
and corrected speed. From this figure we can see 
that increasing inlet pressure is destabilizing as 
is increasing inlet temperature with respect to 
corrected speed but not with respect to physical 
speed. Decreasing incidence was found to be de- 
stabilizing as seen from the two sets of data 
shown. With the second stator closed the inci- 
dence angle will be smaller than that with the 
first stator open. 

An indication of choke flutter was obtained 
while accelerating the engine. As this is done 
the overall compressor pressure ratio increases. 
However, it was observed from the third rotor 
aerodynamic stage data that the stage pressure 
ratio was decreasing as the engine was accelerated 
indicating stage choking. This was observed to 
occur much sooner than any strain gage activity. 

The effect of inlet temperature and inlet 
pressure on the flutter frequency is shown in 
Fig. 14. Here we see that the flutter frequency 
decreases as the pressure increases or as the 
temperature decreases . This is indicative of a 
density effect on the flutter frequency. Assuming 
that choke flutter is an instability such as a 
relaxation limit cycle amenable to definition by 
an equation "of the Van der Pol type it can then 
be assumed that the negative work during an os- 
cillatory-cycle is a function of the air density. 
Also in evidence in Fig. 14 is an indication of a 
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change in aeromechanics! mode of four of the flut- 
ter data points. For these four points the most 
responsive blade was different from the other ten 
data points. The phenomena of mode shape changes 
has also been observed in the stall flutter data 
on one occasion but at that time it was attributed 
to an unintentional bleed from a video view port. 
No explanation has been found for the change in 
the turbojet data. 


Turbojet Program - Stall Flutter 

Prior to the turbojet aeromechanical program 
an aerodynamic program was conducted using the same 
model of engine. In this program the stators were 
strain gaged and monitored. During a part of this 
program while acquiring research data a fourth 
stator vane instability was encountered unexpect- 
edly. While taking data at 105.5 percent of cor- 
rected design speed at inlet conditions of 138 kPa 
(20 psia) and -29° C (-20° F) , some action was ob- 
served on the fourth stators that was due to either 
buffeting or flow separation. At 106 percent cor- 
rected speed an instability developed in the sec- 
ond bending mode. Since there was no previous doc- 
umentation on this instability its investigation 
was included in the aeromechanical program. Ten 
data points were obtained, all above 104 percent 
corrected design speed and at inlet conditions from 
138 kPa to 200 kPa (20 to 29 psia) and -29° to 2° C 
(-20° to 35° F). 

Only two strain gages were functional for all 
of the data points obtained, thus the strain gage 
data is somewhat limited. For all but one data 
point the flutter frequency was the same on both 
strain gages although the maximum amplitudes did 
vary. The flutter frequencies were nonengine 
order and varied from 9.47 to 9.91 engine orders. 
The effect of inlet pressure and temperatureson 
the flutter boundaries is shown in Figs. 15 and 
16. As can be seen from these, increasing inlet 
pressure is, as usual, destabilizing. However, 
increasing inlet temperature was found to be stabi- 
lizing which is the inverse of the trend observed 
on the rotating stages. This type of trend is 
indicative of a density effect. The effect of 
inlet pressure and temperature on the flutter fre- 
quency is seen in Fig. 17. Here we see the usual 
trend due to increasing inlet pressure, a reduction 
in the flutter frequency. Increasing inlet tem- 
perature tended to also reduce the flutter fre- 
quency In much the same way as observed in the 
fan stall flutter data. Analysis of the stage 
aerodynamic data did not indicate choking, thus 
the instability is probably a form of stall flut- 
ter. Since this stator has a negative incidence 
angle at the tip the phenomena observed was prob- 
ably negative incidence stall flutter. 

Turbojet Engine — System Mode Instability 

One of the principle priorities of the turbo- 
jet aeromechanical program was to thoroughly 
document a second rotor system mode instability 
that was known to exist on the engine model being 
tested. Unfortunately, the particular engine 
build used was not very sensitive to this insta- 
bility and only one low stress level transient 
data point was obtained. However, some insight 
into this type of instability was obtained. 


The system mode instability was first ob- 
served during the engine's development program and 
was resolved by closing the first stator to de- 
crease the incidence on the second rotor. This 
instability is of interest because it is not com- 
monly seen and thus not well understood. It 
occurred in the second system mode which involves 
a coupling of above shroud first bending and first 
torsion. Stress levels of the order of ±100 MPa 
(±14.5 ksi) were observed in the development pro- 
gram. The instability has two branches as seen 
in Fig. 18. In this program the lower branch was 
not found. The upper branch was found and at the 
time unknowingly transcended at 96 percent of 
physical design speed (82.7% corrected speed). 

It was anticipated that some rather severe 
stress levels would occur once the instability 
boundary was penetrated. This did not occur. The 
instability was approached by doing a slow ac- 
celeration to 100 percent design speed with nominal 
inlet conditions of 200 kPa (29 psia) and 121° C 
(250° F) . Near 96 percent physical speed a slight 
action on the strain gages was noted, less than 
±30 MPa (±4.3 ksi). The excursion was slowed in- 
anticipation of an instability. However, the 
strain gages quickly subsided and nothing more was 
noted for the remainder of the excursion. The 
recorded strain gage data was analyzed after the 
test run and it became apparent only then that the 
slight strain gage activity noted had indeed been 
an instability that had been transversed. 

Analysis of the strain gage data indicated a 
nonintegral order response over a narrow speed 
range. The maximum stress level observed on any 
strain gage was ±24.8 MPa (±3.6 ksi) at 1370 Hz. 
This frequency is not a fundamental blade fre- 
quency but is associated with a second rotor sys- 
tem mode. An expanded spectrum plot of the maxi- 
mum strain gage signal is shown in Fig. 19, which 
shows that the principal response of the strain 
gage is quite removed from an engine order re- 
sponse. A plot showing stress level and fre- 
quency versus physical engine speed is shown in 
Fig. 20. Prior to penetrating the instability 
there was only a slight nonintegral order re- 
sponse at 1374 Hz. At 95.9 percent physical speed 
the stress level began to rise and reached a maxi- 
mum at 96 percent physical speed at a frequency of 
1370 Hz. Shortly thereafter, the stress levels 
began to decrease and the frequency increase. At 
97 percent physical speed the instability was 
cleared. The decrease in frequency noted in this 
data is probably due to increased damping as the 
stress amplitudes increased. This is simply des- 
cribed by the equation: 

m f = io n (l - C 2 ) (5) 

The damping coefficient ? is sensitive to am- 
plitude variations. The two most active strain 
gages were analyzed to determine their phase re- 
lationships. The phase angle between the two at 
the time of maximum stress was 3.07° which meant 
that they were nearly in phase. 

To induce the instability it was necessary 
to open the first stator to increase the incidence 
angle of the second rotor and to also close the 
second stator to increase the stage pressure 
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ratio. In addition the variable geometry schedule 
was shifted to track, within limits, on the stall 
side of the schedule. Thus, the instability was 
induced much like a stall instability. It should 
be noted that there was no evidence of flow sepa- 
ration in any of the data. 

The test conditions used to obtain this data 
were the extreme allowable maximum conditions thus 
it was not possible to further explore this in- 
stability without exceeding established safe en- 
gine practices. 

Concluding Remarks 

As was seen from the data presented several 
types of instabilities were investigated. It was 
found that each type had a unique set of character- 
istics that distinguished it from the others. The 
types of instabilities included two forms of sub- 
sonic stall flutter, choke flutter and a system 
mode instability. A comprehensive compilation of 
aerodynamic and aeromechanical data were obtained 
relevant to each type of instability. This data 
has been and is yet being analyzed by both the 
NASA and the various engine manufacturers involved 
in these programs. 

The Joint Engine System. Research aeromech- 
anical programs are continuing with a number of 
programs planned for the future. Additional data 
will be collected relative to several types of 
aeroelastic instabilities. In addition, the scope 
of the aeromechanical programs has been expanded 
to acquire structural dynamic data relative to a 
wide variety of inlet distortion stimuli. 
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Figure 1. - Compressor performance and stability map. 
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Figure 6. - Schematic of aerodynamic instrumentation locations at each measuring station (viewed looking upstream). 
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Figure 7. - Expanded spectrum plot of fan stall flutter 
at 172 kPa (25 psia) and 171° C (340° F). 
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Figure E. - Expanded spectrum plot of fan stall flutter 
at 172 kPa (25 psia) and 13^ C (275P F). 
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Figure 13. - Choke flutter boundaries as a function of inlet 
temperature and engine speed. 
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Figure 14. - Effect of inlet pressure and temper- 
ature on third rotor choke flutter. 
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Figure 15. - Compressor performance map showing inlet 
pressure and temperature effects on fourth stator in- 
stability boundaries. 
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Figure 16. - Fourth stator instability boundaries as a function of 
inlet pressure, inlet temperature, physical speed and cor- 
rected speed. 
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Figure 17. - Effect of inlet temperature and pressure 
on fourth stator flutter frequency. 
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Figure 18. - System mode instability boundaries as a function of 
inlet temperature and engine speed. 
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Figure 20. - System mode instability stress and 
frequency as a function of engine speed. 
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